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Abstract 

We develop a framework to calculate the density matrix of a pair 
of photons emitted in a decay cascade with partial "which path" am- 
biguity. We describe an appropriate entanglement distillation scheme 
which works also for certain random cascades. The qualitative fea- 
tures of the distilled entanglement are presented in a two dimensional 
"phase diagram". The theory is applied to the quantum tomography 
of the decay cascade of a biexciton in a semiconductor quantum dot. 
Agreement with experiment is obtained. 

1 Introduction 

Two photon cascades with multiple decay paths are candidate sources of 
entangled pairs of photons. Practical implementations of quantum infor- 
mation theory (TJ [2] prefer to deal with qubits that are based only on the 
photons' states of polarization [31 HI [5] . Unfortunately, unless the cascade 
obeys restrictive symmetry conditions, the 2-qubit state associated with the 
polarization of the photon pair is mixed and has negligible entanglement. As 
these symmetry conditions are very hard to achieve, a distillation procedure 
is needed in order to obtain entangled polarization qubits. In this paper 
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we discuss a novel distillation method which proceeds by spectrally filtering 
the photons. The method was successfully implemented by [B] in obtaining 
entangled polarization qubits from the biexciton cascades in semiconductor 
quantum dots [7] . Our aim is to describe a theory that allows one to compute 
the polarization density matrix resulting from a general decay cascade, with 
and without distillation. 

The two photon cascades discussed in this paper are illustrated in Fig. [TJ 
Each of the two decay paths in the figure emits a pair of photons with char- 
acteristic polarization and color. In Fig. [T^i the two decay channels are dis- 
tinguished only by their polarization: One channel gives two horizontally 
polarized photons and the second channel gives two vertically polarized pho- 
tons. In Fig. [Tb the two decay channels are also distinguished by the fre- 
quencies (colors) of the emitted photons. When the difference between the 
photon's frequencies is not too large (compared with the radiative width of 
the photons) we call this "partial which path ambiguity" (as the colors of 
the photons are not a perfect indicator of the decay path). In Fig. [lb the 
outgoing photons are spectrally filtered so that only a fraction of the pho- 
tons, those that do not distinguish between the decay channels, are collected. 
These photons are the ones that have equal probabilities to be emitted in 
either channel. 

The two photons state, emitted by any one of the cascades in Fig. [H is 
a pure entangled state. It is entangled because the quantum decay proceeds 
simultaneously along the two decay channels. This, however, does not imply 
that the associated pair of qubits, describing the state of polarization, are 
entangled. The state of the qubits is obtained from the quantum state of 
the photon field by tracing out all the degrees of freedom of the two photons 
(e.g. colors) save the polarization [S]. This state is in general, mixed, and 
possibly unentangled in contrast with the two photons state which is pure 
and entangled. In fact, partial path ambiguity caused by a detuning A that is 
large compared with radiative life times — normally the smallest energy scale 
in the problem — gives negligible entanglement of the two qubits. Fortunately, 
in this case, the entanglement can be distilled by erasing the "which path" 
information [21 [TU] as indicated in Fig. [TJ(c). In fact, by choosing a sufficiently 
narrow window, one can distill maximally entangled pairs. The price one pays 
is that the probability of finding close to maximally entangled pairs is then 
very small. 

Decay cascades with "partial which-path ambiguity" are naturally found 
in the biexciton radiative cascade of semiconductor quantum dots [TH [T2] . 
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Figure 1: A decay cascade where the excited state, \u), decays to the ground 
state \d), along two decay paths each emitting two photons. The left (right) 
branch emits two photons that are vertically, y (horizontally, x) polarized. 
In a) the intermediate level is degenerate and the cascade has perfect "which 
path ambiguity". In b) the degeneracy of the intermediate level is slightly 
broken by the detuning A = \E X — E y \. This cascade has "partial which path 
ambiguity" . c) shows how the entanglement can be distilled through spectral 
filtering by a window of width w that erases the "which path" information 



In these solid-state devices, there is an additional complication in that the 
energy levels of the cascade, are (correlated) random quantities that undergo 
slow (on the radiative time scale) fluctuations [13]. These arise from ran- 
dom variations in the electrostatic potential in the sample. The ensemble of 
photons emitted by the cascade is then a mixed state. Entanglement may or 
may not not be distilled in the case of general random cascades with large 
fluctuations. However, as we shall see, for a standard model of the random 
biexciton cascade, distillation works even when fluctuations are large [6]. 

Our theory allows one to compute the 4x4 density matrix, p, of the 
two photon polarization from the spectral properties of the cascade. More 
precisely, we shall see that p, is determined by the quantum energies and life 
times of the energy levels in the cascade, their distribution, and the spectral 
width of the filter. All these quantities can either be measured or determined 
by the experimentalist. The theory avoids modelling the radiating system 
and we do not need to write a Hamiltonian for the radiating system. What 
we do need, instead, is a "universal" form for the photon state generated by 
a radiating (dipole) cascade. This 2-photon quantum state depends param- 
eterically on the energies and lifetimes of the cascade. The theory applies 
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irrespective of the nature of the source, be it a quantum dot, an atom, a 
molecule etc. It allows us to calculate the measure of entanglement [TO] for 
a given cascade, with and without distillation. It also allows us to optimize 
the flux of entangled pairs. 

The paper is organized as follows: In section[2]we describe the polarization 
density matrix for cascades with two decay channels. In section [3] we describe 
the state of the emitted photons from the radiative cascade in the dipole 
approximation. We describe the entanglement distillation in section H] and 
the magnitude of the non-diagonal elements. In section [5] we discuss the 
phases of the non-diagonal elements of the density matrix. In section [6] we 
extend the theory to random cascades relevant to biexciton in quantum dots 
and in section [7] we compare our theory with the experimental results of 
Akopian et al. [6]. 

2 The polarization density matrix 

Consider a radiating system, say a quantum dot, inside a micro-cavity fol- 
lowed by an appropriate optical setup for photon collection so that the out- 
going radiation propagates along the positive z axis. The polarization of the 
outgoing photons then lies in the xy plane. The initial state of the system 
at time zero is an excited dot while the photon field is in its vacuum state, 
\u) cg> |0), see Fig. [TJ For times much longer than the decay time of the dot, 
1/r, the dot is in the bottom state and the photon field has a pair of freely 
propagating photons, and the quantum state of the dot and photon filed is 
\d) ® Each decay path emits a pair of photons with a characteristic 
polarization: vertical polarization for the left path and horizontal for the 
right path [T4]. The state of the freely propagating pair of photons is then 
necessarily of the form 

\m) = E / ^ 1 ^ 2 A i o ;i (A; 1 ,A; 2 )e^l^l + ^l)*4 lj .4 2i .|0> . (1) 
j=x,y 

Xj are the branching ratios for the two decay modes, a k ■ is a photon creation 

operator with wave vector k and polarization j. Since a kl -aj. 2 ■ is symmetric 
in ki and k 2 only the symmetric part of the functions /C2) contributes 

to the integral reflecting the fact that photons are Bosons. We denote by a s 
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the symmetrization of a ,i.e. 

a s { kl M) = a{klM)Jra{k2M ' (2) 



Since the initial state was normalized and the evolution is unitary, so is 
the final state 

3 

We are interested in the correlations between the polarizations of two 
photons. This is fully described by the reduced polarization density matrix 
whose entries are given by [15] 



k\,k2 

With given by Eq. [TJ one finds for p 

( |A,| 2 7 \ 





V 7 |A y | 2 / 



7 = K\^ s x \a s y ) (5) 



in the basis \xx), \xy), \yx) and \yy) (x and y denote the state of polarization). 
This special form expresses the fact that the amplitude for all processes 
involving the polarization states \xy) and \yx) vanish. Note that the matrix 
has normalized trace, I7I < | and that the state is mixed for | (a x \a y ) \ < 1. 

The two qubits are maximally entangled when | A^- 1 2 = I7I = |. When 
7 = the polarization state is separable and may be thought of as a classical 
random source of correlated qubits. I7I is a measure of the entanglement 
known as the negativity [Ml [12], (being the negative eigenvalue of the partial 
transposition of p.) 

In the following sections we describe a theory that allows us to compute 
7 as a function of the spectral properties of the cascade. 



3 Photons in the dipole approximation 

To make progress we need to know the functions ctj of Eq. (CQ). For this we 
need to make some assumptions about the nature of the radiating system. 
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Consider sources that are small compared with the wavelength of the radia- 
tion they emit. For such sources the dipole approximation applies. We shall 
further assume that the interaction between the source and radiation field is 
weak so that the rotating wave approximation applies [18]. In this setting, 
which applies to a wide varieties of radiating systems, the function aj can 
be calculated explicitly. For a radiative cascade with a single branch this 
function is given e.g. in [T81 [TJJ]. The case of two branches is then simply a 
weighted superposition, as in Eq. (JT]). 

For each branch the function atj can be expressed in terms of the spectral 
properties of the cascade: Zi = Eg — iTg, £ — x,y, u. Ei is the energy of the 
£-th state (we chose the ground state to have zero energy, E d = 0) and Yg is 
its width0- For a dipole at the origin one has [T9"] : 

ajih, k 2 ) = A(k 2 , Zj) A(h + fc 2 , Z u ) (6) 

where 

A(A,Z) = ^| ) Z = E-zT (7) 

and we use units where h = c = 1. This reduces the problem of computing 
the entanglement 7 of Eq. (JSJ) to computing integrals. 



3.1 The limit of small radiative width 

In most applications, the radiative widths Tg are the smallest energy scale in 
the problem. This is the case for the biexciton decay in quantum dot where 

~ 0.8fieV), the detuning A = \E X — E y \ ~ 27fieV) and the energies of 
the emitted photons are much larger [2U|, Eg — Eg ~ 1.32 eV). 

The smallness of T leads to simplifications in many of the integrals which 
can then be evaluated analytically. For example, A(k, Z) is concentrated 
near E with a width T, so, in the limit that T is small, one makes only a 
small error by replacing \k\ by k. It then follows that, to leading order in 
Y/E 

{ A\A) = J«\A\>»lJ (k _» +p =l (8, 

In general, as in the case of biexciton decay, the two photons emitted in 
each cascade have different colors, namely, 

r«|(K-^-)-(^-^)| (9) 

1 The common convention [19] replaces our T by T/2. 
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This distinguishes the two photons which may therefore be treated as distinct 
particles and one may forget about the symmetrization, Eq. Mathemat- 
ically, this follows from the observation that in computing overlaps, products 
of the form 

A(h, Zj)A*(k 2 , Z^Afa + k 2 , Z u )\ 2 (10) 

are small and can be neglected. 

This allows us to immediately show that the entanglement in a cascade 
with partial which path ambiguity is negligible when A ^> T. This follows 
from 

(of K 5 > » (a x \a y ) « (^p) 2 « ^ (11) 

The different colors of the emitted photons resolve the "which path ambigu- 
ity". This mixes the two qubits and essentially kills the entanglement. 

4 Entanglement distillation 

The entanglement can be distilled by selecting those photons which does not 
betray the decay path [6]. Let us denote the average intermediate states 
(exciton) energy as 2E = E x + E y . The first emitted photon do not betray 
the decay path provided one only looks at energies near E u — E. Similarly, 
the second photon does not betray the decay path provided one only collects 
photons with energies E. 

In practice, the distillation is done by filtering the photons through a 
spectral window function. The photons are detected only if their energy is 
either within a window of width w centered at E u — E or within one centered 
at E. This is implemented by a monochromator (an energy filter) which 
transmits only a selected part of the emission spectrum [6]. 

Most photon pairs, are of course, lost in the distillation process. Roughly, 
the fraction of photon pairs that are filtered is of the order 0(w/A), as 
most photons lie in the window of width A + 0(T). One might worry that, 
to be effective, the window must be of the order of the radiative life-time, 
w = O(r). If that was the case, only a very small fraction of the photon 
pairs could be distilled. As we shall see, however, this is not the case. In fact, 
one may choose w = 0(A) so a substantial fraction of the photons will be 
distilled while obtaining considerable entanglement. The price one pays for 
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filtering is that the source is not "on demand" but rather a random source 
of entangled photons [5|. 

The filtering process is represented in the theory by a projection operator 
W, which eliminates from a Fock state all photons spare those whose energy 
lies within appropriate energy windows (irrespective of polarization). Here we 
are only interested in the two photon component of the state after filtration. 
Therefore, one can effectively express the action of the operator W as 

W : aj(k u k 2 ) —> w{k l )w{k 2 )a j {ki, k 2 ) (12) 

where w(k) is the step function 

f 1 \k-{E u -E)\<w/2 
w(k) = I 1 \k-E\<w/2 (13) 

y otherwise 

Evidently, If is a projection operator, i.e. W 2 = W. The identity W = 1 
(w = oo) represents no filtering. 

The distillation succeeds with probability pw = (ip\ W and produces 
the (normalized) filtered state 

\*<) - ^= ("J 
\/Pw 

The filtered, or distilled, density matrix can be computed from the distilled 
state. In particular, for the entanglement, as measured by 7 of Eq. (jHJ) we 
find 

Id = K\— — — , Pw = J2 l^'l ( a J W 4 ( 15 ) 
Pw FT,, 

This reduces the problem to computing integrals, where we account for W 
by summing only the appropriate wavevectors. Fig. [2] shows the probability 
to detect a pair of photons and 7^ of the distilled state, as function of the 
width of the spectral window w. To plot the figure we use parameter values 
corresponding to biexciton decay in a quantum dot. As one expects, the 
entanglement is a decreasing function of w, (for a window of zero width one 
gets a maximally entangled state). On the other hand, the probability that 
the detection succeeds is, of course, an increasing function of the width. 

The qualitative behavior of entanglement distillation can be gleaned by 
inspection of Fig. [3J The function a x is concentrated in a small neighborhood 
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Figure 2: The probability that distillation succeeds is an increasing function 
of the window w shown in the red dotted curve. The entanglement of the 
distilled pair is a decreasing function of w shown in the blue curve. The pair 
is maximally entangled at w — 0. The plot is drawn for parameter values 
corresponding to biexciton decay where A/T « 17. 

of size 0(T) near the point of intersection of the green and blue curve. Sim- 
ilarly, a y is concentrated near the intersection of the purple and blue curve. 
For example, the fact that the entanglement without distillation is small, 
Eq. (llljl . follows from the little overlap of (a x \a y ) each of which is concen- 
trated near a different point. Due to distillation, only amplitudes contained 
in the intersection of red squares are collected. This does two things. It 
decreases the numerator in Eq. f|T5|) which is bad. However, it also decreases 
the denominator which is good. This decrease is much more significant and 
consequently the entanglement increases. 

Perhaps the most interesting things one learns form Fig. [3] is how wide 
does a window have to be to betray the "which path information". This 
happens when the window contains the points of intersections, either red 
with purple, or red with green. If the window does not contain these points 
the path is not betrayed. Since the points have a small size, 0(T), this 
implies that the size of the optimal window is of the scale of the detuning, 
w = A — O(r). Because this window is not small, the probability that the 
distillation succeeds is not very small either. 
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Figure 3: The wave function a x is large near the intersection of the diagonal 
strip, expressing total energy conservation, with the line which represents 
the x decay path where the second photons has energy E x . Similarly, the 
function a y is large near the intersection of the diagonal strip with the line 
representing the y decay path where the second photon has energy E y . The 
cross represents a filter of narrow width. The filter collects photons that are 
contained in the intersection of the cross. 
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5 The phase problem 



From the perspective of quantum information theory the phases in the density 
matrix are gauge dependent quantities (as they are not invariant under local 
unitary operations (21]). Even if Alice and Bob fix the projectors |0) (0| and 
|1) (1|, there is still a freedom to choose the phases of the states 

| a) ® 1 6) -> e^e'*" \a) ® \b) , a,be{0, 1}. (16) 

We refer to this as gauge freedom. Such a transformation will change the 
phases of the non-diagonal entries of the density matrix 

p ab , cd e^ + ^-^p ab , cdl a, 6, c, d e {0, 1} (17) 

in the computational basis. Any reasonable measure of entanglement, and 
in particular | of Eq. (jSJ), is clearly independent of the choice of gauge. 

Quantum tomography [22J is an algorithm to convert 16 measurements 
to the 16 (complex) entries of the density matrix p (describing the ensemble) 
[23] . This means that any quantum tomography algorithm must fix both the 
projectors representing the "computational" basis and fix the gauge. 

In the context of photon polarization the canonical choice (which we used 
throughout this paper) of the "computational" basis is the projectors asso- 
ciated with the x and y linear polarizations. The remaining gauge freedom 
is 

a kj -> e^a kj , a{- -> e'^a^ j = x,y (18) 
for two orthogonal polarizations j. Fixing the right circular polarization by 

«U = < 19 > 

fixes the gauge since 

a kR e l ^a kR , a\ R e-**a\ R , (20) 

requires that all the 0's are the same. This then fixes the phase of 7. 

The phase of 7, which was measured in [6], [2H [25] have, so far, not been 
explained by a theoretical model. In the following, we calculate this phase 
and describe the physical information that is encoded in it. 

Eq. ( flBl determines 7 in terms of the product of the branching ratios 
X X X* and the overlap (af | W \oty)- In the next section, we shall show that 
for a decay cascade with partial which path ambiguity and time-reversal 
invariance, A^A* > 0. It then follows that the phase of 7 is fully determined 
by the phase of (af | W \ cty). 
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5.1 Branching ratios 

The branching amplitudes Xj of Eq. ([!]) are proportional to the appropriate 
dipole matrix elements 

X x = {(u\X\j){j\X\d) 

X y = ((u\Y\j)(j\Y\d) (21) 

where ( is a overall normalization constant and X is the x component 
of the (possibly multi-electron) position operator and similarly Y is the y 
component of the position operator. It follows that 

X* x X y = |C| 2 (d\ X \x) (x\ X \u) (u\ Y \y) (y\ Y \d) (22) 

Observe first, that this quantity is independent of the gauge choice of the 
states \£) of the source, as every ket is paired with the corresponding bra. 
We shall now show that in the case that all the states \i) are non-degenerate, 
there is a choice of gauge so that each matrix element in the product is real. 

Let T denote the antiunitary operator associated with time reversal [261 
[27], i.e. 

(T£\Tk) = (k\£) (23) 

In the case that all states \£) are non-degenerate T \£) = e l/3e \£). By changing 
the gauge to \£) — > e J/3f//2 \£), one sees that \£) may be chosen so that T\£) = 
\£). The position operator is evidently even under time reversal e.g. T*XT = 
X. Plugging this in the definition of the dipole matrix elements we see that 

(£\ X \k) = (£\ T*XT \k) = (T£\ X \Tk) = (k\ X \£) (24) 

We have therefore shown that A* X y is a real number. We shall now show 
that under rather weak continuity assumptions, it must actually be positive. 
X* x X y is a function of the spectral properties of the cascade, and in particular, 
it is a function of the detuning A. It has the same sign as A varies so long 
as the two decay paths are indeed effective (none of the branching ratios, Xj, 
vanishes). It is therefore enough to determine the sign at a single point. We 
shall now give a symmetry argument that at A = one has X x X y > 0. 

Assume that the degeneracy A = is a consequence of (possibly approx- 
imate) rotational symmetry in the x-y plane of the radiating system, (this 
is the case in quantum dots). Since the sign of the product of dipole ma- 
trix elements changes continuously as the Hamiltonian is deformed, we may 
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compute the sign for the case where the rotational symmetry is exact. In 
this case, as the initial, non-degenerate, state \u) must be a state of angular 
momentum about the z-axis. Since angular momentum is conserved the 
final two photon state must also be a state of zero angular momentum about 
the z-axis. 

In this case, perfect which path ambiguity and zero angular momentum 
imply that the state of the outgoing photons is 

\m) = J ^ ate, fc 2 y (|feil+|fe2l)t (<X,. + <A„r) i°) > (as) 

By comparing this with Eq. [T] one easily sees that this state implies Xj = 
l/v2- Hence, j — ~ in Eq. (jHj), which determines the sign of the product 

KXy > 0. 



5.2 The role of the complex pole 

It follows from the analysis above that the phase of 7 is the same as the phase 
of (a x \ W \ a y ). The latter is determined by a two-dimensional integration of 
the function 



w(k 1 )w(k 2 ) A(k\ + k 2 , Z u ) A(k 2 , Zj)A(k 2 , Zj> 



(26) 



The first three factors are positive, and weigh the integrand. The third factor 
may be interpreted as guaranteeing approximate conservation of total energy 
since 



lim 

r^o 



A{k,Z) =S(\k\-E) 



(27) 



This means that to leading order in T the matrix elements of p are determined 
by a one-dimensional integral over k of the function 



w(E u - k)w(k) A(k, Zj)A(k, Zj 



(28) 



The phase of 7 is governed by the phase of A(k, Z x )A*(k, Z y ) which is rep- 
resented graphically in Fig. 4 



Strong filtering: 

Suppose the filtering window W is very narrow with width r < w <C A. 
The window restricts the domain of integration to a very narrow region. The 
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detection probability is small and scales linearly with the window's width 
Pw = 0(j£). The phase of 7 is n—AT/A and the magnitude is approximately 



The state is close to a maximally entangled state. This gives the upper left 
triangle of Fig. [5] and the left part of Fig. [2j 

No filtering: 

No filtering corresponds to W — 1 and a width w = 00. Exact degeneracy, 
A = 0, gives a maximally entangled state, 7 = 1/2. (The two arrows in 
Fig. 4 are complex conjugates.) 

When A > T the integrals are dominated by the neighborhood of the 
poles at k = E x ,E y . The off-diagonal element is almost purely imaginary 
and 7 = («■£)• This accounts for the lower right hand triangle of Fig. \5\ 
and the right hand part of Fig. [21 

6 Random cascades 

Radiative cascades with partial which path ambiguity are found naturally 
in semiconductor quantum dots where \u) is the ground state of a bound 
state of a pair of two electrons and two holes (a biexciton). The states \x) 
and \y) are the ground and first excited state of the bright exciton (a bound 
electron- hole pair). The state \d) describes an empty quantum dot. In this 
case the energies and states slowly fluctuate due to electrostatic transients 
in the semiconductor hosting the quantum dots. In typical cases [6], the 
fluctuations are large (comparable to A) and slow (the time scale is much 
longer than 1/r). 

One can model the situation by letting the spectral properties of the 
cascades, namely Z^s), be appropriate functions of a random variables s, 
with measure dP(s). The two-photon state of Eq. (TjQ), \ip(s)) is then a random 
variable. 

The photon field \4>(s)) of Eq. depends on s through the fluctuating 
complex energies, Z e , of Eq. (JOJ) - The 2 photon state emitted by the dot is 
then described not by a pure state but rather by a density matrix 




(29) 




(30) 
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(a) 



(b) 



(c) 




Figure 4: The phase is determined by the product of the two complex num- 
bers (k — Z x ){k — Z*). The numbers are represented as arrows pointing from 
k to the location of the complex energies Z E (blue arrow for Z x and red ar- 
row for Zy). The location of k is restricted by the window function W. In 
figure (a) the levels are detuned and the spectral window is smaller than the 
detuning. In figure (b) the levels are degenerate and in figure (c) the levels 
are detuned and the window is larger than the detuning. 
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5 10 15 20 25 30 35 

spectral window (in units of T) 



Figure 5: The phase of 7 as a function of the (centered) normalized spec- 
tral window width w/T and the detuning A/T for parameters appropriate 
to biexciton in a quantum dot. The triangular reagin above the diagonal 
represents the situation of a narrow filter where 7 ~ —1/2. The triangular 
region below the diagonal is where the window is large, 7 is small and es- 
sentially purely imaginary. At the bottom of the figure the detuning is small 
and 7 ~ 1/2. 



The probability to distill a state describing a specific random event is 

p(s) = (ij(s)\Wms)). (31) 
Therefore, the probability to distill a photon pair is given by 

p(W\ Pr ) = J p(s)P(s)ds (32) 

Similarly, the value of the distilled 7^ is given by averaging over s 

^ fdP(s) P (s) (33) 



where 7(s) is given by substituting \ip(s)} in Eq. f fT5l ). 

In general, large spectral fluctuations can potentially destroy the distil- 
lation based on fixed spectral windows. This, for example, is the situation 
when the values of E u , E x and E y are independent random variables, as il- 
lustrated in Fig. [6]^a). In the figure, the amplitudes which betray the which 
path ambiguity penetrate the filter. When the fluctuations are smaller than 
A one can remedy this by choosing a sufficiently small spectral window. This 
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Figure 6: The gray shapes represent the fluctuations of the cascade. The 
lines diagonal represent conservation of total energy and correlate the two 
gray shapes. They fluctuate as well. The rectangle represents the filtering. 
In a) some of the events penetrate the filter and thus betray the which path 
ambiguity. In b) the fluctuations do not betray the which path ambiguity 
and entanglement can be filtered. 
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situation is illustrated in Fig. 0(a). However, when the fluctuations are larger 
than A distillation becomes impossible. 

A scenario which allows for distillation also when the fluctuations are 
larger than A is illustrated in Fig. 0(b). 

At first, one may think that the second "good" scenario is contrived and 
would not naturally occur. In fact, this is not the case and this scenario is 
the one that describes biexciton drift. The reason is that the energies E u , 
E y and E x are not independent random variables, as in Fig. [6](a), but rather 
dependent random variables. This is because for a biexciton, E u = E x +E y —B 
where B is the biexciton binding energy [28J which is typically more than 
two orders of magnitude smaller than the exciton energy and its dependence 
on s can be safely ignored. A model for a fluctuating spectral diagram is 
then 

E x ^E x + s, E y ^E y + s, E u ^E u + 2s (34) 

This indeed leads to a scenario like the one illustrated in Fig[6](b). 

We note that the biexciton drift described above has only little effect on 
the calculation of 7^ described in section \5\ and Fig This results from the 
rapid decrease of the probability of detection p(s) from its maximal value at 
s = 0. This can be seen in Fig. [6] which plots p(s)dP(s) for finite value of V. 




Drift, (in units of F) 

Figure 7: The relative "weight" of the density matrix p(s). The weight is 
given by A(s) = P(s)Tr(W p(s)) . The plot is renormalized to yield A(0) = 1. 
The plot is obtained with the experimental values as in Sec. 



18 



7 Comparison with experiment 



We now turn to compare the theoretical calculation with the experimental 
data as described by Akopian et al. [6] 

The parameters used in the theory were measured independently!!: r x « 
T y = 0.8 ± 0.2\pieV], T u « 2T X and E x ^ E y = 1.2S[eV), E u = 2.55[eV\, 
A = 27±3|/teV]. The window that was used in the experiment was of width 
w = 25 ± lOf/ieV]. 

The distribution P(s) of the spectral shift was evaluated from the mea- 
sured spectral lines. It was rather wide, with full width middle height of 
about 50[/ieV]. With the values listed above the probability of detection 
Tr(Wp(s)) falls much faster then the distribution P(s) as a function of \s\, 
to half its value at \s\ ~ 5r ~ 8[/ieV]. 

The numerically calculated contribution to the filtered state (i.e. prob- 
ability of detection times probability distribution for s) as a function of the 
spectral drift s for the above parameters is displayed in Fig. 

When we come to compare the theory with the experimental results, we 
must take into account the measurement error of the tomography, as well as 
the errors on the parameters A, F,s and w (the QD or model parameters). 
These are displayed in Fig. [HJ The measured phase in the experiment was 
— 110° ±17° where we have taken into account the effect of the beam split- 
ter. The beam splitter induces a phase shift of 180° between the X and Y 
polarizations of the reflected photon only (this was ignored in Akopian et al. 
[6], where the phase was reported 70°). This is compared to the theoretical 
result —160° ± 45°, which shows a reasonable fit with the experiment. 



8 Conclusion 

We described a framework to calculate the density matrix of a pair of pho- 
tons emitted in a decay cascade with partial which path ambiguity, encoded 
in the energies of the emitted photons. We showed that one can distill the 
entanglement by selecting only the photons possessing "which path" ambi- 
guity and discuss how this distillation by spectral filtering affect the phase 
of the non-diagonal elements of the two photon density matrix. We showed 
that spectral filtering is quite robust and protected from fluctuations in the 

2 We are using the half-width at half maximum (HWHM) convention, while in [6l [19] 
the radiative width is given according to the FWHM convention. 
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Figure 8: Comparison between the experimental results of Akopian et al. 
and the theory. The graph shows a theoretical calculation of the phase as 
a function of the window, w and the detuning, A, both in units of T. The 
calculation uses the experimentally measured parameters w, A, T and s . 
The uncertanties in these parameters result in an area (rather than a point) 
indicated by the error bars. The possible theoretical values of the phase are 
in the area which is bounded by these error bars. This values are compared 
to the experimentally measured phase and error, which is represented in the 
color bar to the left. 

level's energies as long as these fluctuations are correlated. Our calculations, 
quantitatively describe measurements performed on semiconductor quantum 
dots. 
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